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Abstract 

Measurements of the jet energy calibration and transverse momentum resolution in 
CMS are presented, performed with a data sample collected in proton-proton colli- 
sions at a centre-of-mass energy of 7 TeV, corresponding to an integrated luminosity 
of 36 pb~ . The transverse momentum balance in dijet and 7/Z+jets events is used to 
measure the jet energy response in the CMS detector, as well as the transverse momen- 
tum resolution. The results are presented for three different methods to reconstruct 
jets: a calorimeter-based approach, the "Jet-Plus-Track" approach, which improves 
the measurement of calorimeter jets by exploiting the associated tracks, and the "Par- 
ticle Flow" approach, which attempts to reconstruct individually each particle in the 
event, prior to the jet clustering, based on information from all relevant subdetectors. 
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1 Introduction 

Jets are the experimental signatures of quarks and gluons produced in high-energy processes 
such as hard scattering of partons in proton-proton collisions. The detailed understanding of 
both the jet energy scale and of the transverse momentum resolution is of crucial importance 
for many physics analyses, and it is an important component of the systematic uncertainty. 
This paper presents studies for the determination of the energy scale and resolution of jets, per- 
formed with the Compact Muon Solenoid (CMS) at the CERN Large Hadron Collider (LHC), 
on proton-proton collisions at ^/s = 7 TeV, using a data sample corresponding to an integrated 
luminosity of 36 pb^ . 

The paper is organized as follows: Section |2] describes briefly the CMS detector, while SectionlS] 
describes the jet reconstruction methods considered here. Sections |4] and |5] present the data 
samples and the experimental techniques used for the various measurements. The jet energy 
calibration scheme is discussed in Section |6] and the jet transverse momentum resolution is 
presented in Section [7| 

2 The CMS Detector 

A detailed description of the CMS detector can be found elsewhere HI. A right-handed coordi- 
nate system is used with the origin at the nominal interaction point (IP). The x-axis points to the 
center of the LHC ring, the y-axis is vertical and points upw^ard, and the z-axis is parallel to the 
counterclockwise beam direction. The azimuthal angle (p is measured w^ith respect to the x-axis 
in the xy-plane and the polar angle 9 is defined with respect to the z-axis, while the pseudora- 
pidity is defined as ?/ = — In [tan (0/2)]. The central feature of the CMS apparatus is a super- 
conducting solenoid, of 6 m internal diameter, that produces a magnetic field of 3.8 T. Within 
the field volume are the silicon pixel and strip tracker and the barrel and endcap calorimeters 
{\tj\ < 3), composed of a crystal electromagnetic calorimeter (ECAL) and a brass/scintillator 
hadronic calorimeter (HCAL). Outside the field volume, in the forward region (3 < |//| < 5), 
there is an iron /quartz-fibre hadronic calorimeter. The steel return yoke outside the solenoid 
is instrumented with gaseous detectors used to identify muons. The CMS experiment collects 
data using a two-level trigger system, the first-level hardw^are trigger (LI) ||2l and the high-level 
software trigger (HLT) ||3J. 

3 Jet Reconstruction 

Jets considered in this paper are reconstructed using the anti-Zcj clustering algorithm IH with 
a size parameter R = 0.5 in the y — (p space. In some cases, jets with a size parameter R = 0.7 
are also considered. The clustering is performed by four-momentum summation. The rapidity 
y and the transverse momentum pj- of a jet w^ith energy E and momentum p = {'Px,'Py,Pz) 



are defined as y = j In f ^^^^ j and pj = J pi + pi respectively. The inputs to the clustering 
algorithm are the four-momentum vectors of detector energy deposits or of particles in the 
Monte Carlo (MC) simulations. Detector jets belong to three types, depending on the w^ay the 
individual contributions from subdetectors are combined: Calorimeter jets, Jet-Plus-Track jets 
and Particle-Flow jets. 

Calorimeter (CALO) jets are reconstructed from energy deposits in the calorimeter towers. 
A calorimeter tower consists of one or more HCAL cells and the geometrically correspond- 
ing ECAL crystals. In the barrel region of the calorimeters, the unweighted sum of one single 
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HC AL cell and 5x5 ECAL crystals form a projective calorimeter tower. The association between 
HCAL cells and ECAL crystals is more complex in the endcap regions. In the forward region, 
a different calorimeter technology is employed, using the Cerenkov light signals collected by 
short and long quartz readout fibers to aid the separation of electromagnetic and hadronic 
signals. A four-momentum is associated to each tower deposit above a certain threshold, as- 
suming zero mass, and taking as a direction the tower position as seen from the interaction 
point. 

Jet-Plus-Track (JPT) jets are reconstructed calorimeter jets whose energy response and res- 
olution are improved by incorporating tracking information, according to the Jet-Plus-Track 
algorithm [5J. Calorimeter jets are first reconstructed as described above, and then charged 
particle tracks are associated with each jet, based on the spatial separation between the jet axis 
and the track momentum vector, measured at the interaction vertex, in the Tj — (p space. The 
associated tracks are projected onto the front surface of the calorimeter and are classified as 
in-cone tracks if they point to within the jet cone around the jet axis on the calorimeter surface. 
The tracks that are bent out of the jet cone because of the CMS magnetic field are classified as 
out-of-cone tracks. The momenta of charged tracks are then used to improve the measurement 
of the energy of the associated calorimeter jet: for in-cone tracks, the expected average energy 
deposition in the calorimeters is subtracted and the momentum of the tracks is added to the jet 
energy. For out-of-cone tracks the momentum is added directly to the jet energy. The Jet-Plus- 
Track algorithm corrects both the energy and the direction of the axis of the original calorimeter 
jet. 

The Particle-Flow (PF) jets are reconstructed by clustering the four-momentum vectors of 
particle-flow candidates. The particle-flow algorithm (HIZl combines the information from all 
relevant CMS sub-detectors to identify and reconstruct all visible particles in the event, namely 
muons, electrons, photons, charged hadrons, and neutral hadrons. Charged hadrons, electrons 
and muons are reconstructed from tracks in the tracker. Photons and neutral hadrons are re- 
constructed from energy clusters separated from the extrapolated positions of tracks in ECAL 
and HCAL, respectively. A neutral particle overlapping w^ith charged particles in the calorime- 
ters is identified as a calorimeter energy excess with respect to the sum of the associated track 
momenta. The energy of photons is directly obtained from the ECAL measurement, corrected 
for zero-suppression effects. The energy of electrons is determined from a combination of the 
track momentum at the main interaction vertex, the corresponding ECAL cluster energy, and 
the energy sum of all bremsstrahlung photons associated with the track. The energy of muons 
is obtained from the corresponding track momentum. The energy of charged hadrons is deter- 
mined from a combination of the track momentum and the corresponding ECAL and HCAL 
energy, corrected for zero-suppression effects, and calibrated for the non-linear response of the 
calorimeters. Finally, the energy of neutral hadrons is obtained from the corresponding cal- 
ibrated ECAL and HCAL energy. The PF jet momentum and spatial resolutions are greatly 
improved with respect to calorimeter jets, as the use of the tracking detectors and of the high 
granularity of ECAL allows resolution and measurement of charged hadrons and photons in- 
side a jet, which together constitute ~85% of the jet energy. 

The Monte Carlo particle jets are reconstructed by clustering the four-momentum vectors of 
all stable (ct > 1 cm) particles generated in the simulation. In particular, there are two types of 
MC particle jets: those where the neutrinos are excluded from the clustering, and those where 
both the neutrinos and the muons are excluded. The former are used for the study of the PF 
and JPT jet response in the simulation, while the latter are used for the study of the CALO 
jet response (because muons are minimum ionizing particles and therefore do not contribute 
appreciably to the CALO jet reconstruction). 



The Particle-Flow missing transverse energy (£ j), which is needed for the absolute jet energy 
response measurement, is reconstructed from the particle-flow candidates and is defined as 

Ut = ~T^{^i sin Oj cos (piX + £, sin 0, sin (pif) = £ jX + £ yf, where the sum refers to all candi- 

dates and x, y are the unit vectors in the direction of the x and y axes. 

4 Event Samples and Selection Criteria 

In this Section, the data samples used for the various measurements are defined. In all samples 
described below, basic common event preselection criteria are applied in order to ensure that 
the triggered events do come from real proton-proton interactions. First, the presence of at least 
one well-reconstructed primary vertex (PV) is required, with at least four tracks considered in 
the vertex fit, and with |z(PV)| < 24 cm, where z(PV) represents the position of the proton- 
proton collision along the beams. In addition, the radial position of the primary vertex, p(FV), 
has to satisfy the condition |0(PV) < 2 cm. 

Jet quality criteria ("Jet ID") have been developed for CALO jets fSl and PF jets f9\, which 
are found to retain the vast majority (> 99%) of genuine jets in the simulation, while reject- 
ing most of the misidentified jets arising from calorimeter and/or readout electronics noise in 
pure noise non-collision data samples: such as cosmic-ray trigger data or data from triggers on 
empty bunches during LHC operation. Jets used in the analysis are required to satisfy proper 
identification criteria. 

4.1 Zero Bias and Minimum Bias Samples 

The zero bias and minimum bias samples are used for the measurement of the energy clustered 
inside a jet due to noise and additional proton-proton collisions in the same bunch crossing 



(pile-up, or PU), as described in Section 6.2 The zero bias sample is collected using a random 
trigger in the presence of a beam crossing. The minimum bias sample is collected by requiring 
coincidental hits in the beam scintillating counter |3| on either side of the CMS detector. 

4.2 Dijet Sample 

The dijet sample is composed of events with at least two reconstructed jets in the final state 
and is used for the measurement of the relative jet energy scale and of the jet pj resolution. 
This sample is collected using dedicated high-level triggers w^hich accept the events based on 
the value of the average uncorrected pj {pr not corrected for the non-uniform response of the 
calorimeter) of the two CALO jets with the highest pj (leading jets) in the event. The selected 
dijet sample covers the average jet pj range from 15 GeV up to around 1 TeV. 

4.3 7+jets Sample 

The 7+jets sample is used for the measurement of the absolute jet energy response and of 
the jet Pt resolution. This sample is collected with single-photon triggers that accept an event 
if at least one reconstructed photon has pr > 15 GeV. Offline, photons are required to have 
transverse momentum pj > 15 GeV and \rj\ < 1.3. The jets used in the 7+jets sample are 
required to lie in the \rj\ < 1.3 region. The 7+jets sample is dominated by dijet background, 
where a jet mimics the photon. To suppress this background, the following additional photon 
isolation and shower-shape requirements [lOJ are applied: 

• HCAL isolation: the energy deposited in the HCAL w^ithin a cone of radius R = 0.4 
in the f] — (p space, around the photon direction, must be smaller than 2.4 GeV or less 
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than 5% of the photon energy (E-y); 

• ECAL isolation: the energy deposited in the ECAL within a cone of radius R = 0.4 
in the t] — (p space, around the photon direction, excluding the energy associated 
with the photon, must be smaller than 3 GeV or less than 5% of the photon energy; 

• Tracker isolation: the number of tracks in a cone of radius R = 0.35 in the f] — (p 
space, around the photon direction, must be less than three, and the total transverse 
momentum of the tracks must be less than 10% of the photon transverse momentum; 

• Shower shape: the photon cluster major and minor must be in the range of 0.15-0.35, 
and 0.15-0.3, respectively. Cluster major and minor are defined as second moments 
of the energy distribution along the direction of the maximum and minimum spread 
of the ECAL cluster in the rj — (p space; 

The selected 7+jets sample covers the p^ range from 15 GeV up to around 400 GeV. 

4.4 Z(;f+;/~)+jets Sample 

The Z(/^(+/^^)+jets sample is used for the measurement of the absolute jet energy response. It 
is collected using single-muon triggers with various pj thresholds. Offline, the events are re- 
quired to have at least two opposite-sign reconstructed global muons with j)j > 15 GeV and 
\}]f'\ < 2.3 and at least one jet with |//| < 1.3. A global muon is reconstructed by a combined 
fit to the muon system hits and tracker hits, seeded by a track found in the muon systems only. 
The reconstructed muons must satisfy identification and isolation requirements, as described 
in Ref. |1TT| . Furthermore, the invariant mass M^^, of the tw^o muons must satisfy the condi- 
tion 70 < Mpifi < 110 GeV. Finally, the reconstructed Z is required to be back-to-back in the 
transverse plane with respect to the jet with the highest pj- \A(p{Z,jet) \ > 2.8rad. 

4.5 Z(e+e~)+jets Sample 

The Z(e+e^)+jets sample is used for the measurement of the absolute jet energy response. It 
is collected using single-electron triggers with various pj thresholds. Offline, the events are 
required to have at least two opposite-sign reconstructed electrons with pj > 20 GeV in the 
fiducial region \f]\ < 1.44 and 1.57 < \r]\ < 2.5 and at least one jet with \ri\ < 1.3. The 
reconstructed electrons must satisfy identification and isolation requirements, as described in 
Ref. [11 J. Furthermore, the invariant mass Mee of the electron-positron pair must satisfy the 
condition 85 < Mee < 100 GeV. Finally, the reconstructed Z is required to be back-to-back in 
the transverse plane with respect to the jet with the highest pj- \A(p{Z,iet) \ > 2.7rad. 

5 Experimental Techniques 

5.1 Dijet ;7T-Balancing 

The dijet pj-balancing method is used for the measurement of the relative jet energy response 
as a function of rj. It is also used for the measurement of the jet pj resolution. The technique 
was introduced at the CERN pp collider (SPPS) [12] and later refined by the Tevatron exper- 
iments IIT3llT4t . The method is based on transverse momentum conservation and utilizes the 
pr-balance in dijet events, back-to-back in azimuth. 

For the measurement of the relative jet energy response, one jet (barrel jet) is required to lie 
in the central region of the detector (|/y| < 1.3) and the other jet (probe jet) at arbitrary ?/. 
The central region is chosen as a reference because of the uniformity of the detector, the small 
variation of the jet energy response, and because it provides the highest jet p^-reach. It is 
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also the easiest region to calibrate in absolute terms, using 7+jet and Z+jet events. The dijet 



calibration sample is collected as described in Section 4.2 Offline, events are required to contain 
at least two jets. The two leading jets in the event must be azimuthally separated by A^ > 
2.7rad, and one of them must lie in the \r]\ < 1.3 region. 

The balance quantity B is defined as: 
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where p^j*^ is the average pr of the two leading jets: 
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The balance is recorded in bins of f^f"'"'"' and pj"^. In order to avoid a trigger bias, each p^'^ bin 
is populated by events satisfying the conditions of the fully efficient trigger with the highest 
threshold. 

The average value of the B distribution, (B), in a given /yP''"^'^ and p'^'^ bin, is used to determine 
the relative response TZrei- 
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The variable T^-j-ei defined above is mathematically equivalent to {pj ) / {pj ) for narrow 
bins of p'^^ . The choice of pfp*" minimizes the resolution-bia s effec t (as opposed to binning in 



pbarrei ^ which leads to maximum bias) as discussed in Section 



5.4.1 



below. 



A slightly modified version of the dijet p^-balance method is applied for the measurement of 
the jet Pj resolution. The use of dijet events for the measurement of the jet pj resolution was 
introduced by the DO experiment at the Tevatron IITSi w^hile a feasibility study at CMS was 
presented using simulated events IIT6| . 

In events with at least two jets, the asymmetry variable A is defined as: 
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w^here pj and pj refer to the randomly ordered transverse momenta of the two leading jets. 
The variance of the asymmetry variable c^ can be formally expressed as: 
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If the two jets lie in the same r] region, pj = {pj ) = {pj ) and cr{pj) = CT{pj ) = cr{pj ). 
The fractional jet pj resolution is calculated to be: 
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^^ = Vi.^. (6) 

Pt 

The fractional jet pj resolution in the above expression is an estimator of the true resolution, in 
the limiting case of no extra jet activity in the event that spoil the pj balance of the two leading 
jets. The distribution of the variable A is recorded in bins of the average p^ of the two leading 

jets, p^'^ = \Pt + Pt ) /^' ^^'^ ^^^ variance is proportional to the relative jet pj resolution, 
as described above. 

5.2 7/Z+jet ;7T-Balancing 

The 7/Z+jet px-balancing method is used for the measurement of the jet energy response and 
the jet Pt resolution with respect to a reference object, which can be a 7 or a Z boson. The 
Pt resolution of the reference object is typically much better than the jet resolution and the 
absolute response Rai,s is expressed as: 

P? 

VT 

The absolute response variable is recorded in bins of p^' . It should be noted that, because 
of the much worse jet pj resolution, compared to the 7 or Z p^ resolution, the method is not 
affected by the resolution bias effect (see Section 5.4.1|, as it happens in the dijet pj-balancing 



method. Also, for the same reason, the absolute response can be defined as above, without the 
need of more complicated observables, such as the balance B or the asymmetry A. 

5.3 Missing Transverse Energy Projection Fraction 

The missing transverse energy projection fraction (MPF) method (extensively used at the Teva- 
tron Iil3il ) is based on the fact that the 7, Z+jets events have no intrinsic £ 7 and that, at parton 
level, the 7 or Z is perfectly balanced by the hadronic recoil in the transverse plane: 

,^ 7/Z I ,T* recoil r, /q\ 

Vt +Pt =0- (8) 

For reconstructed objects, this equation can be re-written as: 

R^,ZPV'^ + RrecoUPV'""' = -It, (9) 

where R^,z and Rjecoii are the detector responses to the 7 or Z and the hadronic recoil, respec- 
tively. 

Solving the two above equations for Rjecoii gives: 
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^recoil = K-y,z H , ^7,T = ^UVT- (10) 
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This equation forms the definition of the MPF response Rmpf- The additional step needed 
is to extract the jet energy response from the measured MPF response. In general, the recoil 
consists of additional jets, beyond the leading one, soft particles and unclustered energy. The 
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relation Ruadjet = ^recoil holds to a good approximation if the particles, that are not clustered 
into the leading jet, have a response similar to the ones inside the jet, or if these particles are in 
a direction perpendicular to the photon axis. Small response differences are irrelevant if most 
of the recoil is clustered into the leading jet. This is ensured by vetoing secondary jets in the 
selected back-to-back 7, Z+jets events. 

The MPF method is less sensitive to various systematic biases compared to the 7, Z p^-balancing 
method and is used in CMS as the main method to measure the jet energy response, while the 
7, Z pj-balancing is used to facilitate a better understanding of various systematic uncertainties 
and to perform cross-checks. 

5.4 Biases 

All the methods based on data are affected by inherent biases related to detector effects (e.g. pj 
resolution) and to the physics properties (e.g. steeply falling jet pj spectrum). In this Section, 
the two most important biases related to the jet energy scale and to the pj resolution measure- 
ments are discussed: the resolution bias and the radiation imbalance. 

5.4.1 Resolution Bias 

The measurement of the jet energy response is alw^ays performed by comparison to a reference 
object. Typically, the object with the best resolution is chosen as a reference object, as in the 
7/Z+jet balancing where the 7 and the Z objects have much better pj resolution than the jets. 
However, in other cases, such as the dijet p^-balancing, the two objects have comparable reso- 
lutions. When such a situation occurs, the measured relative response is biased in favor of the 
object with the worse resolution. This happens because a reconstructed jet pr bin is populated 
not only by jets whose true (particle-level) pr lies in the same bin, but also from jets outside the 
bin, w^hose response has fluctuated high or low. If the jet spectrum is flat, for a given bin the 
numbers of true jets migrating in and out are equal and no bias is observed. In the presence of a 
steeply falling spectrum, the number of incoming jets with lower true pj that fluctuated high is 
larger and the measured response is systematically higher. In the dijet pj-balancing, the effect 
described above affects both jets. In order to reduce the resolution bias, the measurement of the 
relative response is performed in bins of p^p*^, so that if the two jets have the same resolution, 
the bias is cancelled on average. This is true for the resolution measurement w^ith the asymme- 
try method w^here both jets lie in the same r] region. For the relative response measurement, the 
two jets lie in general in different rj regions, and the bias cancellation is only partial. 

5.4.2 Radiation Imbaiance 

The other source of bias is the px-imbalance caused by gluon radiation. In general, the mea- 
sured p^-imbalance is caused by the response difference of the balancing objects, but also from 
any additional objects with significant pj. The effect can by demonstrated as follow^s: an esti- 
mator 'JZ™'^"'^ of the response of an object w^ith respect to a reference object, is T?.""^"^ = pj/pj 
where pj and Pj are the measured transverse momenta of the objects. These are related to 
the true pr iPj^rlf) through the true response: pj = R*™'^ • p'/"^ and pj = R*™/ ■ Pj"g<r. In 
the presence of additional hard objects in the event, p^""^ = Pj"lf — ^Pt, where Apr quantifies 
the imbalance due to radiation. By combining all the above, the estimator TZ""^"^ is expressed 
as: TZ"""' = Rf^"<^/Rtrue fi _ i^^^i^true\ -^Us relation indicates that the pr-ratio between two 
reconstructed objects is a good estimator of the relative response, only in the case where the 

true 



additional objects are soft, such that Ap^/p^"^, — > 0. 
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The above considerations are important for all pj-balancing measurements presented in this 
paper (the dijet pr-balancing and the 7/Z+jet pr-balancing), both for the scale and the reso- 
lution determination. Practically, the measurements are performed with a varying veto on an 
estimator of fl^'"'^ = Apj/p^''^, and then extrapolated linearly to fl^''"'^ = 0. For the dijet pj- 

balancing, the estimator of a*'"'' is the ratio pj /p^^, while for the 7, Z+jet pj-balancing it is 

the ratio p['^^/pj'^. 

6 Jet Energy Calibration 

6.1 Overview of the Calibration Strategy 

The purpose of the jet energy calibration is to relate, on average, the energy measured for the 
detector jet to the energy of the corresponding true particle jet. A true particle jet results from 
the clustering (with the same clustering algorithm applied to detector jets) of all stable particles 
originating from the fragmenting parton, as well as of the particles from the underlying event 
(UE) activity. The correction is applied as a multiplicative factor C to each component of the 
raw^ jet four-momentum vector p™'" (components are indexed by fi in the foUow^ing): 

p';' = c-p;^^. (11) 

The correction factor C is composed of the offset correction Coffset/ the MC calibration factor 
CmC/ arid the residual calibrations Qei and Cabs for the relative and absolute energy scales, 
respectively. The offset correction removes the extra energy due to noise and pile-up, and the 
MC correction removes the bulk of the non-uniformity in 7/ and the non-linearity in p^. Finally, 
the residual corrections account for the small differences between data and simulation. The 
various components are applied in sequence as described by the equation below: 

w^here pj is the transverse momentum of the jet after applying the offset correction and pj is 
the Pt of the jet after all previous corrections. In the following sections, each component of the 
jet energy calibration will be discussed separately. 

6.2 Offset Correction 

The offset correction is the first step in the chain of the factorized corrections. Its purpose 
is to estimate and subtract the energy not associated with the high-pj scattering. The excess 
energy includes contributions from electronics noise and pile-up. In CMS, three approaches 
are followed for the offset correction: the jet area, the average offset and the hybrid jet area 
methods. 

6.2.1 Jet Area Method 

Recent developments in the jet reconstruction algorithms have allow^ed a novel approach for 
the treatment of pile-up HzHlHl: for each event, an average pr-density p per unit area is esti- 
mated, which characterizes the soft jet activity and is a combination of the underlying event, 
the electronics noise, and the pile-up. The two latter components contaminate the hard jet 
energy measurement and need to be corrected for with the offset correction. 
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The key element for this approach is the jet area Aj. A very large number of infinitely soft four- 
momentum vectors (soft enough not to change the properties of the true jets) are artificially 
added in the event and clustered by the jet algorithm together with the true jet components. The 
extent of the region in the y — (p space occupied by the soft particles clustered in each jet defines 
the active jet area. The other important quantity for the pile-up subtraction is the pj density 
p, which is calculated w^ith the kj jet clustering algorithm |[191 - |2T| w^ith a distance parameter 
R = 0.6. The kj algorithm naturally clusters a large number of soft jets in each event, which 
effectively cover the entire y — <p space, and can be used to estimate an average pr-density The 
quantity p is defined on an event-by-event basis as the median of the distribution of the variable 
PjJ Aj, w^here ] runs over all jets in the event, and is not sensitive to the presence of hard jets. 
At the detector level, the measured density p is the convolution of the true particle-level activity 
(underlying event, pile-up) with the detector response to the various particle types. 

Based on the know^ledge of the jet area and the event density p, an event-by-event and jet-by-jet 
pile-up correction factor can be defined: 



Carea(pr.A>P)=l 



[p - {p\J^)) ■ Aj 
y,raw 

Ft 



(13) 



In the formula above, (pue) is the pj-density component due to the UE and electronics noise, 
and is measured in events with exactly one reconstructed primary vertex (no pile-up). FigurejT] 
shows the PF p^-density p, as a function of the leading jet pr in QCD events and for various 
pile-up conditions. The fact that p does not depend on the hard scale of the event confirms that 
it is really a measure of the soft jet activity. Finally, the density p shows linear scaling properties 
with respect to the amount of pile-up. 
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Figure 1: Pile-up and underlying event PF p^-density p, as a function of the leading jet pj 
in the QCD multijet sample for various pile-up conditions (here Npv denotes the number of 
reconstructed vertices, and A denotes the unit area in the y — ^ space). 



6.2.2 Average Offset Method 

The average offset method attempts to measure the average energy due to noise and pile-up, 
clustered inside the jet area, in addition to the energy associated w^ith the jet show^er itself. The 
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measurement of the noise contribution is made in zero bias events by vetoing those that pass 
the minimum bias trigger. In the remaining events, the energy inside a cone of radius R = 0.5 
in the t] — (p space is summed. The measurement is performed in cones centered at a specific 
rj bin and averaged across (p. The noise contribution is found to be less than 250 MeV in pj, 
over the entire f/ range. The total average offset (over the entire dataset) is determined from 
inclusive zero bias events (with no veto on minimum bias triggers) and is classified according 
to the number of reconstructed vertices. Figure [2] shows the average offset pj as a function of 
rj and for different pile-up conditions. The calorimetric offset pr show^s strong variations as a 
function of tj, w^hich follow the non-uniform particle response in the calorimeter, while for PF 
candidates, the offset pr is more uniform versus rj. The higher measured offset pj for the PF- 
candidates is due to the much higher response with respect to the pure calorimetric objects. The 
observed ^/-asymmetry is related to calorimeter instrumental effects. For the highest number of 
vertices, in particular, the asymmetry is also of statistical nature (the adjacent points are highly 
correlated because at a given t] a large fraction of the energy in a cone of R = 0.5 also ends 
up in overlapping cones). Figure |3] shows the breakdown, in terms of PF candidates, of the 
average offset pj in events with one PU interaction, as measured in the data and compared to 
the MC prediction. The slight asymmetry observed in the MC is due to the asymmetric noise 
description in the specific version of the simulation. The average offset in pj scales linearly with 
the number of reconstructed primary vertices, as shown in Fig.H The linear scaling allows the 
expression of the jet offset correction as follows: 

Offset W/Pt ''^FV) - 1 -rPE^ ' (^4) 

Ft 

where 0(?/) is the average pr due to one pile-up event, p™'" is the pj of the uncorrected jet, and 
Npv is the number of reconstructed primary vertices. The average offset method can be applied 
to jet algorithms that produce circular jets, while the quantity 0{r]) scales to larger cone sizes in 
proportion to the jet area. It should be noted that, in both the average offset subtraction and in 
the jet area method, the noise contribution and the UE are not subtracted. Because of the good 
description of the noise contribution in the simulation, the noise is taken into account w^ith the 
MC-based correction. 

6.2.3 Hybrid Jet Area Method 

The measurement of the average offset presented in the previous paragraph confirms the //- 
dependence of the offset energy. This is explained by the fact that the measured offset is the 
convolution of the pile-up activity with the detector response. In order to take into account the 
^/-dependence, a hybrid jet area method is employed: 

Ft 



In Eq. (15 >, the pj density p and the corresponding density due to the UE, (pue) are constants 
over the entire t] range. The multiplicative factor /3(7/) corrects for the non-uniformity of the 
energy response and is calculated from the modulation of the average offset in pj (Fig.[2]|. 

In the case of PF jets, the response variation versus ?/ is relatively small and the hybrid jet 
area method is found to be in excellent agreement with the average offset method. Figure [4] 
shows the average offset in pr as a function of the number of reconstructed primary vertices, 
for the three different methods (jet area, average offset, hybrid jet area). It can be seen that the 
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Figure 2: Average offset in pj, as a function of t], measured in minimum bias events for dif- 
ferent pile-up conditions (categorized according to the number Npv of reconstructed primary 
vertices). Left: CALOjets. Right: PFjets. 
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Figure 3: Breakdown of the average offset pr, in terms of the PF candidates, as a function of 
rj, for events with one PU interaction. Data are show^n by markers and MC is shown as filled 
histograms. 



differences betw^een the jet area method and the average offset method are entirely due to the 
response dependence on t]. The hybrid jet area method is chosen for the pile-up correction of 
PF jets. In the case of CALO jets, and also JPT jets (initially reconstructed as CALO jets), the 
response variation versus tj shows dramatic changes and neither the simple jet area nor the 
hybrid jet area methods are able to reproduce the average offset measurement. Therefore, for 
CALO jets, the average offset method is the one chosen for the pile-up correction. 

6.2.4 Offset Uncertainty 

The uncertainty of the offset correction is quantified using the jet area method. Specifically, the 
quantities p and (pue) in Eq. ([l3| are varied independently and the resulting shifts are added 
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for the jet area, the average offset, and the hybrid jet area methods in 6 different r] regions. In the 
y-axis title, PFAK5 denotes the PF jets reconstructed with the anti-Zcj algorithm with distance 
parameter R = 0.5. 



in quadrature. The event pj-density p uncertainty is estimated as 0.2 GeV per unit jet area 
and per pile-up event. This uncertainty is based on the maximum slope difference between 
the jet area and the average offset methods, and the residual non-closure in the average offset 
method. The UE p^-density (pue) uncertainty is estimated as 0.15 GeV per unit jet area, based 
on the differences observed between the QCD multijet and Z-i-jets samples, and on the effective 
difference when applied in the inclusive jet cross-section measurement. Figure [5] shows the 
uncertainty of the offset correction, as a function of jet p^ and the number of primary vertices. 

6.3 Monte Carlo Calibration 

The MC calibration is based on the simulation and corrects the energy of the reconstructed 
jets such that it is equal on average to the energy of the generated MC particle jets. Simulated 
QCD events are generated with PYTHIA6.4.22 [22J, tune Z2 (the Z2 tune is identical to the 
Zl tune described in [23 1 except that Z2 uses the CTEQ6L PDF, while Zl uses CTEQ5L) and 
processed through the CMS detector simulation, based on GEANT4 p4l . The jet reconstruc- 
tion is identical to the one applied to the data. Each reconstructed jet is spatially matched in 
the y] — <p space with a MC particle jet by requiring AR < 0.25. In each bin of the MC particle 



transverse momentum p ^ , the response variable 7?. = jfj^° / p ^ and the detector jet jfj^° are 
recorded. The average correction in each bin is defined as the inverse of the average response 
Cmc(Pt™) ~ <fe' ^^'^ ^^ expressed as a function of the average detector jet pj < pj™ >. 
Figure |6] shows the MC jet energy correction factor for the three jet types, vs. ?/, for different 
corrected jet pr values. Figure [t] shows the average correction in |?^| < 1.3, as a function of the 
corrected jet pj. 

Calorimeter jets require a large correction factor due to the non-linear response of the CMS 
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Figure 5: Offset jet-energy-correction uncertainty as a function of jet pr. Left: CALO jets. Right: 
PF jets. 



calorimeters. The structures observed at I?/ 1 ~ 1.3 are due to the barrel-endcap boundary and 
to the tracker material budget, which is maximum in this region. The fast drop observed in the 
endcap region 1.3 < |f/ 1 < 3.0 is due to the fact that the jet energy response depends on energy 
rather than on jet pr. For higher values oi \r]\ more energy corresponds to a fixed pj value 
E ^ Pt ■ cosh(f/), which means that the jet response is higher and the required correction factor 
is smaller. The structure observed Sit\r]\ ~ 3.0 coincides with the boundary between the endcap 
and the forw^ard calorimeters. Finally, in the region \rj\ > 4.0, the jet energy response is lower 
because parts of the jets pointing toward this region extend beyond the forward calorimeter 
acceptance. 

The track-based jet types (JPT and PF) require much smaller correction factors because the 
charged component of the jet shower is measured accurately in the CMS tracker w^hich extends 
up to I //I = 2.4. The fast rise of the correction factor for JPT jets in the region 2.0 < \f]\ < 2.5 
is explained by the fact that part of the jets lying in this region extends beyond the tracker 
coverage. For PF jets, the transition beyond the tracker acceptance is smoother because the PF 
candidates, which are input to the clustering of PF jets, are individually calibrated prior to the 
clustering. While both PF jets and JPT jets exploit the tracker measurements, the JPT jets require 
lower correction in the region \ri\ < 2.0 because the tracker inefficiency is explicitly corrected 
for by the JPT algorithm. In the forward region (I?/! > 3.0) all three jet types converge to simple 
calorimetric objects and therefore require almost identical corrections. 

The default MC calibration is derived from the QCD sample and corresponds to a jet flavour 
composition enriched in low-pj- gluon jets. The jet energy response and resolution depend on 
the fragmentation properties of the initial parton: gluons and heavy-flavour quarks tend to 
produce more particles with a softer energy spectrum than light quarks. The investigation of 
the jet energy response of the various flavour types, for the different jet reconstruction tech- 
niques, is done with MC matching between the generated particle jet and the reconstructed jet. 
For each MC particle jet, the corresponding parton is found by spatial matching in the ^ — (p 
space. Figure IS] shows the response of each flavour type (gluon, b-quark, c-quark, uds-quark), 
as predicted by PYTHIA6 (Z2 tune), in the region \f]\ < 1.3, normalized to the average re- 
sponse in the QCD flavour mixture. The QCD flavour composition varies significantly with 
jet Pt, being dominated by gluon jets at low pj and by quark jets at high pr- Calorimeter jets 
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Figure 7: Monte Carlo jet-energy-correction factors for the different jet types, as a function of 
jet Pt- 



show strong dependence on the flavour type with differences up to 10%. This is attributed to 
the non-linear single-particle response in the calorimeters. For the track-based reconstructed 
jets, the flavour dependence is significantly reduced and not larger than 5% and 3% for JPT 
and PF jets respectively. The ability to measure precisely the charged particle momenta in the 
tracker reduces the contribution of calorimetry at low jet pj- In all jet types, the jets originated 
from a light quark (u/d/s) have a systematically higher response than those from the other 
flavours, which is attributed to the harder spectrum of the particles that are produced in the 
fragmentation process. For comparison, Fig.[9]shows the flavour dependent response ratio of a 
different fragmentation model (HERWIG++) with respect to PYTHIA6. 
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6.4 Relative Jet Energy Scale 

6.4.1 Measurement 

The dijet pr-balance technique, described in Section |5| is used to measure the response of a jet 
at any rj relative to the jet energy response in the region \f]\ < 1.3. Figure 10 shows example 



distributions of the balance quantity B for PF jets in two pseudorapidity bins. Figure 11 shows 
the relative response as a function of rj in the range 100 GeV < pf^ < 130 GeV. Ideally, the 
relative response of the corrected jets in the simulation should be equal to unity. However, 



because of the resolution bias effect (Section 5.4.1 1, the relative response in the simulation is 
found to deviate from unity by an amount equal to the resolution bias. The comparison of 
the data w^ith the MC simulations implicitly assumes that the resolution bias in the data is the 
same as in the simulation. This assumption is the dominant systematic uncertainty related to 
the measurement of the relative response with the dijet balance method. 



In order to reduce the radiation bias (Section 5.4.2 1, a selection is applied on the ratio oc 



Jet3 



Pj I p^j^ and the nominal analysis value is a < 0.2. The residual relative correction calculation 
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Figure 10: Example distributions of the dijet balance quantity for PF jets in two // regions. 



is done in three steps: first, the ;/ -symmetric part, Csym, is measured in bins of |?/|, in order to 
maximize the available statistics, with the nominal requirement a < 0.2. Then, a correction 
factor /Cj-ad is applied to take care of the extrapolation to a = 0, and finally the asymmetry in 
?// '^r(|'/|)/ is taken into account. The residual correction for the relative jet energy scale is 
formally expressed below^: 



Crel{±V) 



^raddj/D-Csymd?/!) 



(16) 



The Csym component is defined by comparing the relative response in data and MC simulations 



Qymd'/I) = ( ^a<o.2 
\ data , 



(17) 



Pt 



averaged over the entire pj range. This is justified by the fact that no statistically significant 
p7-dependence is observed in the comparison betw^een data and simulation. 

Since the additional radiation and the UE are not perfectly modeled in the simulation, a correc- 



tion needs to be applied by extrapolating to zero third-jet activity, as discussed in Section 5.4.2 
The radiation correction fc^arf is defined as: 



^rad 




(18) 



Figure [12] (left) shows the radiation correction that needs to be applied to the measurement at 
the working point a. < 0.2. The correction is negligible in the central region while it reaches the 
value of 3% at larger rapidities. 

The asymmetry of the response in rj is quantified through the variable Ar: 
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Figure 11: Relative jet energy response as a function of t], measured with the dijet balance 
method for CALO, JPT and PF jets respectively. 



M\v\) 



R{+\V\) - R{-\t]\) 
R(+|,7|)+K(-|^|)' 



(19) 



where R(+|?/|) (R( — |f/|)) is the relative response measured in the data at the detector part 



lying in the direction of the positive (negative) z-axis. Figure 12 (right) show^s the measured 
asymmetry. It is found to be similar for the different jet types. 



Figure 13 show^s the final residual correction, as a function of r], for all jet types. This correction 
istypicallyoftheorder of 2-3%, with the exception of the region 2.5 < \t]\ < 3.0 where it reaches 
the value of 10%. The region where the larger discrepancy between data and MC simulations is 



observed (Fig. Ill, coincides w^ith the border between the endcap and the forward calorimeters. 
It has also been observed p5| that the single-particle response show^s similar behavior in this 
region. 



Finally, Fig. 14 demonstrates that the derived residual correction establishes an almost perfect 
agreement between data and simulation. 
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Figure 12: Left: correction /Cj-ad of the relative jet eriergy residual due to initial and final state 
radiation. Right: relative jet energy response asymmetry as a function of jet \r]\, for a. < 0.2. 



6.4.2 Uncertainty 

The dominant uncertainty of the relative residual correction is due to the simulation of the jet 
energy resolution, which defines the magnitude of the resolution bias. The estimate of the sys- 
tematic uncertainty is achieved by varying the jet pj resolution according to the comparisons 
between data and MC simulations shown in Section IT] Other sources of uncertainty, such as 
lack of available events, radiation correction and asymmetry in rj are found to be smaller than 



1%. The total uncertainty of the relative jet energy scale is shown in Fig. 15 as a function of 
the jet \rj\ for two characteristic values of jet pr (50GeV, 200 GeV). The CALO jets have sys- 
tematically larger uncertainty, as opposed to PF jets which have the smallest while the JPT jets 
uncertainty lies between the values for the other two jet types. This pattern is consistent with 
the behavior of the jet energy resolution. Also, it is observed that the relative scale uncertainty 
grows toward larger rapidities because of the larger resolution uncertainty. 
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Figure 13: Relative jet energy residual correction as a function of jet rj for CALO, JPT and PF 
jets respectively. The band shows the uncertainty due to statistics, radiation corrections, and 
asymmetry in tj. 
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Figure 14: Relative response ratio between data and MC simulation before and after the resid- 
ual correction. 
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Figure 15: Relative jet energy residual correction uncertainty, as a function of r] for jet j)j 
50 GeV (left) and pr = 200 GeV (right). 
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6.5 Absolute Jet Energy Scale 

6.5.1 Measurement 

The absolute jet energy response is measured in the reference region \t]\ < 1.3 with the MPF 
method using 7/Z+jets events, and the result is verified with the pr-balancing method. The 
7 or the Z are used as reference objects because their energy is accurately measured in ECAL 
(photon, Z — )■ e+e^) or in the tracker and muon detectors (Z — )■ fi^ji^). Figure 16 shows 
example distributions of the MPF and p^-balancing methods for PF jets in the 7+jet sample. 

The actual measurement is performed only for PF jets because of the full consistency between 
the jet and the Bt reconstruction (both use the same PF candidates as inputs). The absolute 
energy scale of the remaining jet types (CALO, PF) is determined by comparison to the corre- 
sponding PF jet after jet-by-jet matching in the tj — (p space. 
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Figure 16: Example response distributions for PF jets from pj-balancing (left) and MPF (right) 
in the 7+jets sample. 



In the selected 7+jets sample, the presence of a barrel jet ( | f/ 1 < 1.3) recoiling against the photon 
candidate in azimuth by A(p > 2.7 is required. To reduce the effect of initial and final state 
gluon radiation that degrades the jet-photon p^-balance, events containing additional jets with 
Pj > a. ■ pj and outside the AR = 0.25 cone around the photon direction are vetoed. The 
pr-balance and MPF response measurements are performed in the same way w^ith data and 
MC samples with different values of the threshold on a and the data/MC ratio is extrapolated 
to a = 0. This procedure allows the separation of the 7-jet intrinsic p^-imbalance from the 
imbalance caused by hard radiation (Section 5.4.2| |. 

Figure [17| (left) shows the data/MC jet-energy-response ratio, relative to the 7 ECAL scale, ex- 
trapolated as a function of the threshold on the second jet pj- In the p7-balancing method, 
the secondary jet effect is more pronounced because it affects directly the transverse momen- 
tum balance between the photon and the leading jet. In the MPF method, the presence of the 
secondary jet(s) affects the measurement to a lesser extent, and mainly through the response 
difference between the leading jet and the secondary softer jet(s). For loose veto values, the ra- 
tio data/MC in both methods is lower than unity, while the agreement improves by tightening 
the veto. Figure 17 (right) shows the data/MC response ratio after the extrapolation to a = 
for both MPF and pj-balancing methods, as a function of pj. The tw^o measurements are statis- 
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Figure 17: Left: dependence of the data/MC ratio of the jet energy response on the second jet pr 
threshold. Right: data/MC ratio of the jet energy response, after extrapolation to zero second 
jet Pt, as a function of pj. Solid squares and solid circles correspond to the pf-balancing and 
the MPF methods, respectively. 



tically uncorrected to a good approximation and the two sets of points are fitted together with 
a constant value. The fit gives data/MC = 0.985 ± 0.001, relative to the 7 ECAL scale, which 
leads to an absolute response residual correction Cabs = 1/0.985 = 1.015 (Eq. ([12]|), constant in 

Pt- 

In addition to the 7+jets sample, the absolute jet energy response is also measured from the 
Z+jets sample. Figure 18 shows two characteristic response distributions in the 30 GeV < p| < 
60 GeV bin, as an example, measured from the Z(/^+^^ )+jets sample w^ith the pj-balancing and 
the MPF methods. The Z+jets samples cover the pj range from 20 GeV to 200 GeV. 
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Figure 18: Left: jet energy response from Z{]A^ji )+jets pr-balancing in the bin 30 < pf < 
60 GeV. Right: jet energy response from Z(f/+/i^ )+jets MPF in the bin 30 < pf < 60 GeV. 



In order to combine the results from the photon+jet and Z+jet samples, the more precise MPF 
method is employed identically in all relevant samples. Figure Il9] shows the data/MC ratio as 
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Figure 19: Ratio of data over MC for the MPF response, as a function of pj.' in the photon+jet 
sample (circles), Z(e+e^)+iet sample (triangles) and Z(/^(+/^^)+jet sample (squares). 



a function of pj' after correcting for the final and initial state radiation differences between 
data and simulation (extrapolation to a = 0). Although the size of the Z+jets data sample is 
smaller than the 7+jets sample, the results from all samples are in good agreement, within the 
corresponding statistical uncertainties. 

6.5.2 Uncertainty Sources 

The uncertainty of the absolute jet energy scale measurement has six components: uncertainty 
in the MPF method for PF jets, photon energy scale, MC extrapolation beyond the reach of 
the available dataset, offset due to noise and pile-up at low-pr (as discussed in Section 6.2.4| |, 
MC residuals (the level of closure of the MC correction in the MC), and the jet-by-jet matching 
residuals for CALO and JPT jets. 

MPF Uncertainty for PF Jets. The MPF method is affected by several small uncertainties that 
mainly contribute at low pj- flavour mapping, parton-to-particle level sensitivity, QCD back- 



ground, secondary jets, and proton fragments. The various contributions are shown in Fig. 20 



The flavour mapping uncertainty accounts for the response difference between jets in the quark- 
rich 7+jets sample used to measure the absolute jet energy scale, and those in the reference, 
gluon-rich QCD multijet sample. This is estimated from the average quark-gluon response 
difference between PYTHIA6 and HERWIG++ (Fig.|9) in the region 30 - 150 GeV. The latter is 
chosen because it is the pj region best constrained by the available data. For PF jets, the flavour 
mapping uncertainty amounts to ~ 0.5%. 

By definition, the MPF response refers to the parton level because the photon is perfectly bal- 
anced in the transverse plane, against the outgoing partons. However, the default jet energy 
response refers to the particle level, which includes the UE and the hadronization effects. The 
parton-to-particle level response interpretation therefore is sensitive to the UE and the out-of- 
cone showering (OOC). The corresponding uncertainty is estimated from the simulation by 
using jets reconstructed with larger size parameter (R = 0.7, more sensitive to UE and OOC) 
and comparing the extrapolation to the zero secondary jet activity with respect to the nominal 
size parameter (R = 0.5). The resulting uncertainty has a weak pj-dependence and is smaller 
than 0.2%. 
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The dominant background for 7+jets events is the QCD dijet production where one leading jet 
fragments into a hard isolated n^ — > 7 + 7. Such events can alter the measured p^-balance 
because the leading neutral n^ carries only a fraction of the initial parton energy. The QCD 
background uncertainty is estimated by repeating the measurement, using a loose and a tight 
photon identification, and is found to be negligible compared to the current statistical precision. 

The MPF response at low pr is sensitive to the undetected energy that leaks outside the forward 
calorimeter acceptance sA \t]\ > 5. This results in an underestimation of the MPF response, 
compared to the true response. The uncertainty due to the undetected energy is taken from the 
simulation and is estimated to be 50% of the difference between the MPF response and the true 
response. 

The secondary jet activity is found to be significantly different betw^een data and MC, and 
it is corrected by extrapolating the data/MC ratio for the MPF and pj-balance methods to 
zero secondary jet activity. The related uncertainty is estimated as half of the radiation bias 
correction applied to the MPF method. 

Photon Energy Scale Uncertainty. The MPF and p^-balancing methods are directly sensitive to 
the uncertainty in the energy of the 7 used as a reference object. The 7 energy scale uncertainty 
is estimated to be ~ 1% based on studies presented elsewhere 



Monte Carlo Extrapolation. The in situ measurement of the absolute jet energy scale is feasible 
only in the pj range w^here 7+jets data are available. For the current dataset this range extends 
to around 300 GeV. However, the jet pj range probed in the entire dataset is generally more 
than three times higher than in the 7+jets sample. In QCD dijet events, jets as high as pj = 
1 TeV are observed. Because of the absence of data for direct response measurement at high pj, 
the calibration relies on the simulation. Based on the data vs. MC comparison in the region of 
available 7+jets data, conclusions can be drawn for the extrapolation of the jet energy correction 
at the highest jet pr- 

The simulation uncertainty for the high-pj jets arises from two main sources: the single-particle 
response (SPR) and the fragmentation modeling. The former is measured directly in data by 
using isolated tracks and comparing the energy deposited in the calorimeters with the mo- 
mentum measured by the tracker. The currently available measurement Il25l indicates that the 
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Figure 21: Left: sensitivity of the jet energy response in |j/| < 1.3 to the single-particle response 
(SPR) uncertainty. Right: dependence of the jet energy response on the fragmentation model. 
Here AK7PF stands for PF jets reconstructed with the anti-A;^ algorithm with size parameter 
R = 0.7. 



data/MC disagreement is less than 3%. The SPR uncertainty is translated to a jet energy re- 
sponse uncertainty by modifying accordingly the simulation. Figure [2l] (left) shows the impact 
of the SPR uncertainty on the response of the different jet types, in the region \rj\ < 1.3. For 
CALO jets, the induced uncertainty is roughly constant vs. pr and approximately equal to 
2%. The track-based algorithms are less affected at low-p^ by the SPR uncertainty because the 
energy is primarily measured by the tracker. However, as the jet pj increases and the track 
momentum measurement becomes less precise compared to the calorimetric measurement, the 
track-based jet types behave like CALO jets. The transition is smooth and is completed at jet 
PT~300GeV. 

The other source of systematic uncertainty is related to the fragmentation properties, which 
include the parton shower and the hadronization simulation. Since jets are composite objects, 
realized as "sprays" of highly collimated particles, and the calorimeter response is non-linear, 
the jet energy response depends on the number and the spectrum of the particles it consists 
of. The sensitivity to the fragmentation modelling is studied by generating QCD events from 
various MC generators which are then processed by the full simulation of the CMS detector. 
The MC generators employed are: PYTHIA6 (tunes D6T f23] and Z2) and Herwig++||2Z1. Fig- 
ure 21 (right) shows the response ratio of the various models with respect to PYTHIA6,with Z2 
tune, which is the default. The differences between the models are negligible at pj ~ 80 GeV, 
while they grow^ up to 1.5% at low and high jet pj. 

The combined MC uncertainty of the absolute jet energy response due to SPR and fragmenta- 



tion is shown in Fig. 22 



The particle-flow^ algorithm reconstructs individual particles, prior to jet clustering. This al- 
lows the detailed study of the PF jet composition in terms of charged hadrons, photons and 
neutral hadrons. In particular, the jet energy response is closely related to the energy fraction 
carried by the three major composition species. The purpose of this study is to demonstrate 
that the MC simulation is able to describe accurately the PF jet composition observed in data 
and therefore can be trusted to predict the PF jet response in the kinematic regions w^here the 
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Figure 22: Uncertainty of the absolute jet energy response in the region \rj\ < 1.3 related to the 
MC extrapolation for CALO, JPT and PF jets respectively. 



in situ measurement is not possible. 

Figure |23] (left) shows the fraction of jet energy carried by the various particle types. Charged 
hadrons, photons, and neutral hadrons carry ~ 65%, 20%, and 15% of the jet energy respec- 
tively at low jet Pt, as expected from the general properties of the fragmentation process. 
As the jet pj increases, charged hadrons become more energetic and more coUimated, w^hile 
the tracking efficiency and momentum resolution worsen. This increases the probability for a 
charged hadron to leave detectable energy only in the calorimeters and to be classified either 
as a neutral electromagnetic object (photon) or as a neutral hadron. Therefore, for higher jet pj, 
the energy fraction carried by photons and neutral hadrons is increased. The excellent agree- 



ment between data and simulation quantified in Fig. 23 (right) proves that the simulation can 
be safely trusted to predict the absolute jet energy response. 

Jet-by-Jet Matching. Once the jet energy scale is established for PF jets, the estimated un- 
certainties are transfered to the other jet types. This is done by direct jet-by-jet comparison 
between different jet types in the QCD dijet sample. The PF and CALO (JPT) jets are spatially 
matched in the t], (p space by requiring AR < 0.25. For the matched jet pairs the relative re- 
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Figure 23: PF jet composition. Left: energy fraction carried by charged hadrons (CHF), photons 
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sponse of CALO (JPT) jets pj^^'~^ /pj^ (pj /pj^) is measured as a function of pj^ (the study is 
described in detail in Ref. [9|). A cross-check of the direct jet matching is done with a tag-and- 
probe method in dijet events, with the PF jet being the tag object and the CALO/JPT jets being 



the probe objects. The results are summarized in Fig. 24 where the response ratio data/MC of 



the CALO and JPT response relative to the PF jets is shown. The observed disagreement is at 
the level of 0.5%, indicating that the precision of the CALO and JPT calibration is comparable 
to that of the PF jets. The observed 0.5% level of data/MC disagreement is taken into account 
as an additional systematic uncertainty for CALO and JPT jets. 
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6.5.3 Uncertainty 

As described in the previous sections, the absolute jet energy response is measured in situ for 
PF jets with the MPF method in 7/Z+jets events. The systematic uncertainties related to the 



measurement itself are summarized in Fig. 20 The estimation of the systematic uncertainty 



in the kinematic region beyond the reach of the 7+jets sample is based on the simulation and 
its sensitivity to the single-particle response and the fragmentation models. In addition, the 
uncertainty on the 7 energy scale needs to be taken into account since the jet energy response 
is measured relative to the 7 scale. The direct jet-by-jet spatial matching, allows the transfer 
of the PF jet-energy- scale uncertainty to the other jet types (CALO, JPT). Finally, a flavour 
uncertainty is assigned from the response differences betw^een the quark and gluon originated 
jets. These are taken from Fig. |9] and cover the absolute scale uncertainty in physics samples 
with a different flavour mixture than the reference QCD multijet sample. 



Figure 25 shows the absolute energy scale uncertainties for the three jet types, combined with 
the offset correction uncertainty corresponding to the average number of pile-up events in the 
datasets considered for this paper. The low jet pr threshold indicates the minimum recom- 
mended Pt for each jet type: 30 GeV, 20 GeV, and 10 GeV for CALO, JPT, and PF jets respec- 
tively. At low jet Pt the offset uncertainty dominates w^ith significant contribution from the MC 
truth and jet-by-jet matching residuals. At the intermediate jet pj, where enough data for the 
in situ measurements are available, the 7 energy scale uncertainty dominates. At high jet pr, 
the uncertainty due to the MC extrapolation is dominant. Overall, the absolute jet energy scale 
uncertainty for all jet types is smaller than 2% for p^ > 40 GeV. 
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Figure 25: Absolute jet energy scale uncertainty as a function of jet pr for CALO, JPT and PF 
jets respectively. 
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6.6 Combined Jet Energy Correction 

In this section, the combined MC and residual calibration is presented along with the total jet 
energy scale systematic uncertainty. Following Eq. p2] |, the residual corrections for the relative 
and absolute response are multiplied with the generator-level MC correction, while the cor- 



responding uncertainties are added in quadrature. Figure 26 shows the combined calibration 
factor as a function of jet-/^ for p^ = 50, 200 GeV. Because of the smallness of the residual cor- 
rections, the combined correction has the shape of the MC component, shown in Fig. |6| The 



total correction as a function of jet pj is shown in Fig. 27 for various rj values. Figure 128 1 shows 



the total jet energy scale uncertainty as a function of jet pj- At low jet pj the relative energy 
scale uncertainty makes a significant contribution to the total uncertainty w^hile it becomes neg- 
ligible at high Pt- In the forward region, the relative scale uncertainty remains significant in 
the entire pj-range. In general PF jets have the smallest systematic uncertainty while CALO 
jets have the largest. 
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Figure 26: Total jet-energy-correction factor, as a function of jet rj for pj 
Pt = 200 GeV (right). The bands indicate the corresponding uncertainty. 
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Figure 27: Total jet-energy-correction factor, as a function of jet pr for various r] values. The 
bands indicate the corresponding uncertainty. 
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Figure 28: Total jet-energy-scale uncertainty, as a function of jet pj for various rj values. 
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7 Jet Transverse Momentum Resolutions 

In the following sections, results on jet pj resolutions are presented, extracted from generator- 
level MC information, and measured from the collider data. Unless stated otherwise, CALO, 
PF and JPT jets are corrected for the jet energy scale, as described in the previous section. 

The jet pr resolution is measured from two different samples, in both data and MC samples, 
using methods described in Section |5] 

• The dijet asymmetry method, applied to the dijet sample, 

• The photon-plus-jet balance method, applied to the 7+jet sample. 

The dijet asymmetry method exploits momentum conservation in the transverse plane of the 
dijet system and is based (almost) exclusively on the measured kinematics of the dijet events. 
This measurement uses two ways of describing the jet resolution distributions in data and 
simulated events. The first method makes use of a truncated RMS to characterize the core of the 
distributions. The second method employs functional fitting of the full jet resolution function, 
and is currently limited to a Gaussian approximation for the jet pr probability density. 

The 7+jet balance method exploits the balance in the transverse plane between the photon and 
the recoiling jet, and it uses the photon as a reference object whose pr is accurately measured 
in ECAL. The width of the pr/pj distribution provides information on the jet pj resolution 
in a given pj bin. The resolution is determined independently for both data and simulated 
events. The results extracted from 7+jet pj balancing provide useful input for validating the 
CMS detector simulation, and serve as an independent and complementary cross-check of the 
results obtained w^ith the dijet asymmetry method. 

In the studies presented in this paper, the resolution broadening from extra radiation activity 
is removed by extrapolating to the ideal case of a tw^o-body process, both in data and in MC. In 
addition, the data/MC resolution ratio is derived. 

7.1 Monte Carlo Resolutions 

The jet pj resolution derived from generator-level MC information information in the simula- 
tion, serves as a benchmark for the measurements of the jet resolution in collision data samples, 
using the methods introduced above and discussed in the following sections. 

The measurement of the jet pj resolution in the simulation is performed using PYTHIA QCD 
dijet events. The MC particle jets are matched geometrically to the reconstructed jets (CALO, 
JPT, or PF) by requiring their distance int] — (p space to be AR < AR^ax- 

The jet pj response is defined as the ratio pj'^" / pj" where pj'^° and p^" refer to the transverse 
momenta of the reconstructed jet and its matched reference MC particle jet respectively. 

The w^idth of the jet pj response distribution, in a given |//| and p^" bin, is interpreted as the 



generator-level MC jet pr resolution. Figure 29 shows an example of p^'^°/p^'^" distribution for 
CALO jets in If/ 1 < 0.5 and with 250 < pf " < 320 GeV. 

7.2 Dijet Measurements 

The principles of the dijet asymmetry method for the measurement of the jet pj resolution were 
presented in SectionlS] Here, the results of the measurement are presented. 

The idealized topology of two jets w^ith exactly compensating transverse momenta is spoiled 
in realistic collision events by the presence of extra activity, e.g. from additional soft radiation 
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Figure 29: Distribution of the simulated CALO jet response, Pj'^° /p^", in a particular |//| and 
p^"^ range. Fit examples with a Gaussian and a double-sided Crystal-Ball function are shown. 

or from the UE. The resulting asymmetry distributions are broadened and the jet pr resolution 
is systematically underestimated. Other effects can also cause jet imbalance. For example, 
fragmentation effects cause some energy to be showered outside the jet cone ("out of cone 
radiation"). The width of the asymmetry distribution is thus a convolution of these different 
contributions: 



^A — ^intrinsic tp (^imhalance 



(20) 



To account for soft radiation in dijet events, the measurement of the asymmetry in each rj and 
pip"^ bin is carried out multiple times, for decreasing amounts of extra activity, and the jet pj 
resolution is extracted by extrapolating the extra event activity to zero, as discussed in Sec- 
The ratio of the transverse momentum of the third jet in the event over the dijet 



tion 



5.4.2 



Ft 



Jets 



average pj, Pj"""" = p'j' /p!p"^, is used as a measure of the extra activity. The extrapolation 
procedure is illustrated in Fig. |30] (left) for the 120 < p^p*" < 147 GeV bin of PF jets and for the 
corresponding bin of MC particle jets (right). The width of each asymmetry distribution aj^, as 
well as the resolutions obtained using generator-level MC information, are derived based on 
the RMS of the corresponding distributions. Some characteristic example distributions for the 



raw asymmetry are shown for PF jets in Fig. 31 



To account for the particle-level imbalance contribution to the measured jet pj resolution, the 
asymmetry method is applied to the generated MC particle jets. Then the extrapolated particle- 



level resolution is subtracted in quadrature from the measurement. Figure 32 illustrates the 
different steps of the asymmetry procedure for CALO, JPT, and PF jets respectively. The total 
pr resolution derived from the extrapolation of the reconstructed asymmetry is shown in green 
circle, the estimation of the particle-level imbalance resolution from the application to MC par- 
ticle jets is shown in magenta diamond, and the quadrature subtraction to the final asymmetry 
result is shown in blue square. All three can be described by a fit to a variation of the standard 
formula for calorimeter-based resolutions. 
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Figure 30: Examples of extrapolations of v2a_/[ as a function of p^*^ / p^'^ to zero for PF jets (R = 
0.5) in |;y| < 0.5 and 120 < pip"^ < 147 GeV (left). Example of a corresponding extrapolation for 
MC particle jets (right). 



Pt 



= Jsgn(N)-(^) +S--fr-''+C^ 



(21) 



where, N refers to the "noise", S to the "stochastic", and C to the "constant" term. The addi- 
tional parameter M is introduced, and the negative sign of the noise term is allowed, to improve 
the fits to the jet pj resolution vs. pr, for jets that include tracking information (JPT, PF), while 
retaining a similar functional form as the one used for CALO jets. The resolution estimated 
from generator-level MC information is shown in red triangles, and good agreement w^ith the 
result of the asymmetry method is observed. The ratio Maarmmetr^T ^^ obtained as a func- 
tion of Pt, for each jet type and in each rj-bin and is later applied to the data measurement as a 
bias correction. 

Several sources of systematic uncertainties are identified: 

The linear extrapolation at half-the-distance betw^een the standard working point (at pj / p^p"^ = 
0.15) and zero is evaluated, and the difference from the full extrapolation to zero is assigned 
as an uncertainty. The size of the particle-level imbalance is varied by 25% and the impact 
of the measurement is studied when subtracting 75% and 125% of the original particle jet pj 
resolution in quadrature. 

Performing the analysis on simulated events, w^e observe deviations (biases) from the obtained 
and expected values, referred to as "MC closure residuals". A conservative 50% of the MC 

1 -J 1 MC(eenerator— level)— MC(asymmetry) . , , jj-i- i i i- i ,.■ 

closure residuals — — MCfas mmetr ) ^^ taken as an additional relative systematic 

uncertainty, corresponding to the bias correction. By comparing the asymmetry measured in 
data with the expectation from MC simulations, an additional constant term is fitted, describing 
the observed discrepancy betw^een data and simulation, as described below. The statistical 
uncertainty from the fit of the constant term is assigned as a systematic uncertainty. Figure 33 
shows the size of the different systematic uncertainties as a function of pf^ and for a central 
rj bin, for the three jets types. The particle-level imbalance uncertainty is shown in opaque 
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Figure 31: Examples of PF jet asymmetry distributions for \ri\ < 0.5 and a low 

a medium-p^^*^ bin (top right) and a high-p ^""^ bin (bottom), determined from QCD simulation 

(blue histograms) and compared with the result from data (black dots). 



orange, the solid yellow contribution corresponds to the uncertainty from the soft radiation 
variation, and the dashed-red line depicts the impact from the remaining differences in the 
MC closure. The relative uncertainty due to particle-level imbalance is larger for JPT and PF 
jets than for CALO jets because the absolute values of the raw resolutions are significantly 
smaller for JPT and PF, and thus more sensitive to the imbalance subtraction, than in the CALO 
jet case. The dashed blue line shows the contribution of the uncertainty on the additional 
constant term. The total systematic uncertainty for each resolution measurement is obtained 
by summing all individual components in quadrature, and is represented by the grey filled 



area in Fig. 33 The sensitivity of the method to the presence of additional collisions due to 
pile-up has been assessed by applying the measurement to the subsample of the data where 
exactly one primary vertex candidate is reconstructed, and no significant deviations from the 
inclusive measurement are observed. 

The presented measurements of the jet pj resolution, obtained by applying the asymmetry 
method to data, yield systematically poorer resolution compared to the simulation. This dis- 
crepancy is quantified by taking the fits to the MC asymmetry results, fixing all parameters, and 
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Figure 32: Application of the asymmetry method to simulated CALO (top left), JPT (top right), 
and PF jets (bottom) in |?/| < 0.5. The reconstruction-level (green circle) and particle-level 
(magenta diamond) results are shown together with the final measurement (blue square), com- 
pared to the generator-level MC (denoted as MC-truth) derived resolution (red triangle). 

adding in quadrature an additional constant term, as the only free parameter in a subsequent 
fit to the data asymmetry. The fitted additional constant term provides a good characterization 
of the discrepancy, which was verified by several closure tests based on MC. A likely source of 
the discrepancy is an imperfect intercalibration of the CMS calorimeters, w^hich affects analyses 
based on the corresponding datasets. 



The final results are presented in Figs. 34 (for all three types of jets, in the central region) and 
35 (for PF jets in all remaining t] bins). In each case, the solid red line depicts the resolution 
from generator-level MC, corrected for the measured discrepancy between data and simula- 
tion (constant term), and represents the best estimate of the jet pj resolution in data. Conse- 
quently, it is central to the total systematic uncertainty band, drawn in yellow. The uncorrected 
generator-level MC resolution is shown as a red-dashed line for reference. The black dots 
are the bias-corrected data measurements, which are found to be in good agreement with the 
discrepancy-corrected generator-level MC, within the statistical and systematic uncertainties. 
Note in particular that the agreement with the uncorrected generator-level MC resolution is 
considerably worse. 
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Figure 33: Relative systematic uncertainty of the asymmetry method to simulated CALO (top 
left), JPT (top right), and PF jets (bottom) for |?/| < 0.5. 

The dijet data are also investigated within the framework of the unbinned likelihood fit to the 
jet Pj resolution parameterization. This approach is developed in order to provide a cross- 
check of the results. It also serves as a tool for the determination of the full jet pj resolution 
function, once larger collider data samples become available. This method directly takes into 
account biases in the event selection caused by the jet pj resolution and the steeply falling jet 
Pt spectrum. 

At the present stage, the jet pj probability densities are approximated by a truncated Gaussian, 
providing direct correspondence with the binned fits discussed above. The resulting determi- 
nation of the widths of the jet pr resolution (as function of pj and rj) is also affected by the 
soft-radiation and hadronization (out-of-cone) effects. The fitted resolution values are thus ex- 
trapolated to zero-radiation activity. The MC particle-level imbalance is subtracted in quadra- 
ture to correct for effects of hadronization. The method is applied to both data and MC, and the 
results are consistent w^ith the previously discussed binned fits to the asymmetry distributions. 
Namely, poorer resolutions are observed in data compared to the simulation. 
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Figure 34: Bias-corrected data measurements, compared to the generator-level MC (denoted as 
MC-truth) Pt resolution before (red-dashed line) and after correction for the measured discrep- 
ancy between data and simulation (red-solid line) for CALO (top left), JPT (top right), and PF 
jets (bottom) in |//| < 0.5. 
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Figure 35: Bias-corrected data measurements, compared to the generator-level MC (denoted 
as MC-truth) resolution before (red-dashed line) and after correction for the measured discrep- 
ancy between data and simulation (red-solid line), compared to data, for PF jets in different rj 
ranges. 
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7.3 7 + Jet Measurements 

As for dijets, the measurement of the jet pj resolution using the 7+jet pj-balancing, involves 
an extrapolation of the event topology to the ideal case of zero secondary hadronic activity, as 



described in detail in Section 5.4.2 To measure the jet pj resolution from data, the observable 
^{pj /pj) is expressed as: 

^totalKPr ' V t) ~ ^intrinskyV J I Vt ) ® '^imbalance kPt 'Vt)' (^^) 

where the first term (JintrinskiVT ^ P^^) i^ ^^e intrinsic (generator-level MC) resolution of inter- 
est. The second term (Timbalanceip^" / Pj) i^ ^^^ "imbalance" term, arising from the presence of 
secondary jets in an event and from other effects, such as hadronization. The effect of extra jet 
activity is studied as a function of pj I p\. 

The jet pj resolution is measured using two methods. The "direct" method measures the pj 
resolution separately for data and MC, while the "ratio" method is specialized for the data/MC 
ratio. In the direct method, the intrinsic resolution is taken to be independent of pj I pj, while 
the width of the imbalance is assumed to be a first-order polynomial in the pj / pj fraction. 
The intercept "q" gives the limit of zero secondary jet activity and the parameter "m" describes 
the soft-radiation effects: 

^intrinsic\P J /Pj) ^ ''' {■^•^) 

0-imbalanceip?^ / PI) =q + m- pf^ / p}, (24) 

^totalip^T^/pl) = ^J{c' + q^ + lqmp^flp\ + m\p^flp\Y) . (25) 

In order to correct for the particle-level imbalance, the measured pj resolution, using recon- 
structed jets and photons, is fitted with the functional form in Eq. p5] l, while keeping the pa- 
rameter c\ fixed to the value obtained from the fit to the imbalance in MC, using the functional 
form in Eq. (|24]|. A non-zero value of the intercept c\ is expected, w^hich represents the irre- 
ducible imbalance due to, e.g. hadronization and photon resolution effects. The parameter c is 
the measured intrinsic resolution from the reconstructed quantities. The same functional form 
is used in fits to both MC and data samples. 



Figure 36 (left) shows p^j /pj distributions, for different bins of Pj / Pj, in MC simulation. The 
effect of tightening the secondary jet activity is clearly visible as a narrowing of the spread of the 
distributions. The intrinsic resolution on the other hand, is independent of any other activity in 
the event. The evolution of the RMS of the p^-balance distribution in bins of pj /pjis show^n 



in Fig. 36 (right) as the red circular points, for both data and MC samples. The blue square 
points show the intrinsic resolution of PF jets in MC simulation, measured as the RMS of the 
Pt /Pt"- ^^^ grey thicker line shows the total expected resolution defined as the quadrature 
sum of intrinsic (blue square points) and imbalance component (green triangle points), as a 
function of pj /pj. The dotted-red line is the result of the fit, applied to MC simulation, and 
the thinner solid-red line is the functional fit to the data measurement. 

The 7+jet sample has significant contamination from QCD dijet background with a jet misiden- 
tified as a photon. To pass the photon isolation and cluster shape requirements, such a jet must 
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Figure 36: Response distributions in different bins of j)j /pj, and for 140 < Pj < 180 GeV 
(left), components of the jet pr resolution, as a function of pj /pj for PF jets in data and MC 
samples (right). 

be composed of a leading k^ (with n^ — )■ 77) accompanied by a very low hadronic activity. 
These misidentified "photons" have energy-scale similar to the genuine prompt photons, a 
good energy resolution, and can therefore serve as valid reference objects for this analysis. In 
the selected photon sample, the presence of a jet is required, with \ri\ < 1.1, recoiling against 
the photon candidate in azimuth within Acp > 2n/3. 

Distributions of the pj/p^ variable for PF jets are shown in Fig. 37 for data and MC. The dis- 
tributions are not centered around the response of 1.0, because of the impact of soft radiation. 



as illustrated in Fig. 36 



The jet pj resolution obtained from the 7+jet analysis, provides a cross-check of the dijet asym- 
metry results, and a reasonable agreement is observed betw^een the two measurements, as 



shown in Fig. 38 At the current level of statistics, the 7+jet method yields poorer resolutions 
for pt > 150 GeV. 

The complementary "ratio" method is based on taking the ratio of the data and MC intrinsic 
resolutions versus pj I pj before the extrapolation. The intrinsic resolutions are first derived 
in data and MC, by subtracting in quadrature, from the total measured resolutions, the im- 
balance predicted in the simulation. The strength of the method is that the extrapolation fit 
is performed only once, and that the fitted observable, as estimator of the ratio of data and 
MC intrinsic resolutions, is expected to be a constant function of Pj /pj. The intrinsic res- 
olution derived in data is confirmed to be flat vs. pj /pj, as expected, providing a test of 
the procedure. Any deviation from a flat dependence would have indicated a limitation of the 
simulation to model the imbalance in the data. The results for the ratio are consistent with the 
constant pj / pj behaviour and therefore a simple constant fit has been performed in bins of 
Pt and rj. Systematic uncertainties due to variation of the extrapolation fit range and the un- 
certainty of jet energy corrections have been evaluated, adding up to ±(3 — 4)% for the ratio. 
An additional uncertainty of ±(2 — 4)% is assigned to the ratio, due to the assumption that 
the MC simulation correctly models the imbalance in data. This uncertainty is evaluated by 
varying the relevant features at generator level, namely, the modeling of hadronization, treat- 
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Figure 37: Distributions of pr/Pj- in data and MC sample for PF jets in two representative pj 
bins. 

ment of multiple-parton final states and modeling of the /cT-kick. The direct and ratio methods 
are consistent with each other within uncertainties. The statistical errors on the results from 
the ratio method are smaller than from the direct method, because of the fact that in the ratio 
method the imbalance is fixed to the MC -based result, while in the direct method the parameter 
m, describing the imbalance part of the resolution, is free in the fits to data and MC samples. 



The results from the ratio method are compared in Fig. 39 (left) to the results from the dijet 
samples (using the unbinned likelihood fits), as a function of pr, in the barrel region. Since 
the tw^o samples have different jet flavour compositions, the generator-level MC resolutions 
from the corresponding MC samples have been compared and verified that the difference is 
within 3%. The dependence of the data/MC ratio on the flavour difference is expected to be 
~ 1% under a conservative assumption of 30% uncertainty on the modeling of the flavour 
composition in the simulation. This component is included in the comparison of the data/MC 
ratio between the dijet and 7+jet samples. 

Therefore, the two results are compared directly, without applying additional jet composition 
corrections. Both sets of points are fully consistent and have been combined in fits, using a 
constant function of pj for each rj range. The ?/ dependence of the measured intrinsic resolution 



ratios is also shown in Fig. 39 (right) 



The results from the ratio method, as well as the results from the dijet samples (using the un- 
binned likelihood fits), integrated over pj and for all four different \f]\ regions are illustrated 
in Table [T] As before, good agreement is observed between the two methods, within the statis- 
tical and systematic uncertainties. The quoted uncertainties for the dijet analysis combine the 
systematic uncertainties assigned to the extrapolation procedure, the particle-level imbalance 
correction, the jet-energy-scale correction, and the particle-level dijet pj spectrum. 

7.4 Measurement of Jet Resolution Tails 

One of the most promising signatures of physics beyond the standard model involves events 
with multiple jets and a large missing transverse energy £ 7. A background to this signal is 
expected from QCD multijet production where £ j can originate, e.g. from fluctuations in the 
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Figure 38: Jet pj resolutions from 7-l-jet (red triangular points) and dijet asymmetry (blue cir- 
cular points) measurements, for CALO (top left), JPT (top right), and PF jets (bottom). 

detector response to jets. One way to estimate the QCD background in the high ^ 7 signal 
region is to smear particle-level multijet events with parameterizations of the full jet pj resolu- 
tion functions that model both the Gaussian core and the tails of the distributions. It is therefore 
important to quantify the non-Gaussian component of the jet pj resolution in order to predict 
accurately the QCD background. 

Two complementary studies of resolution tails are presented, using dijet and 7-l-jet events. 
For these studies, the focus is on the PF jet reconstruction, since it provides the best jet pj 
resolution and is adopted in the primary physics analyses most sensitive to the impact of the 
jet Pt resolution tails. 

7.4.1 Dijet Asymmetry Measurement 

The full resolution functions can be derived using the generator-level MC information in the 
simulation. To validate the MC simulation description of the pj-resolution tails with the cur- 
rently available data samples, the fractional number of events in the tail regions of the dijet pj 
asymmetry distributions is compared between data and simulation. 

As show^n before (Table ll||, the central-core widths of the response and asymmetry distribu- 
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Figure 39: The ratio of jet pj resolutions in data and MC samples vs. pj, in |?7| < 1.1, from dijet 
and 7+jet samples; a combination fit to both data sets is also shown (left). Results from the 
combination fits vs. pj- in various rj ranges (right). 

Table 1: Ratios of the resolution measured in data and simulation for different jet t] ranges 
using the unbinned likelihood fits in dijet samples; the last column gives the ratio for PF jets 
from 7+jet samples for comparison. Stated are the statistical uncertainties from the fit and the 
upper and lower systematic uncertainties. 



I ?/ 1 bin Ratio C ALO Jets Ratio JPT Jets 



Ratio PF Jets Ratio PF Jets in 7+jet 



0.0-1.1 1.088 ± 0.007+°;°^^ 1.087 ±0.0061°;°^° 1.066 ± 0.007+°;°^^ 

1.1-1.7 1.139 ± 0.019+^:^^^ 1.213 ±0.015lg:{^^J 1.191 ± 0.019+g;g^^ 

1.7-2.3 1.082 ± 0.030+^;^^^ 1.018 ±0.02ll°;°^| 1.096 ± 0.030+g;g| 

2.3-5.0 1.065 ± 0.042+° 1^ 1.068 ±0.036+°; j^^ 1.166 ± 0.050+°;j^^ 



1.07±0.020+°:°^4 
1.10 ± 0.031 +g;g^^ 

1.18 ± 0.062+°;°^^ 



tions differ between data and MC samples. The adopted strategy is to adjust the MC response 
distributions to have the same core resolutions in MC simulations as in data. Then, the fraction 
of events in a given asymmetry w^indow in the tail of the distribution is calculated with both 
data and MC samples. Different tail regions have been studied; here the results for the window 
2.5c7" — 00 are presented. These fractions are observed to depend on the threshold on the third-jet 
Pt, and are therefore extrapolated to zero. The measured ratio between data and MC fractions 
from asymmetry is used to correct the fraction from generator-level MC in the form of a scal- 
ing factor. Since the dijet asymmetry distribution is symmetric by construction, the measured 
tail-scaling factors average over the low and the high response tails. To validate the method 
and to quantify biases caused by the event selection, the extrapolated fractions in the MC sim- 
ulation have been compared to the expectation from the asymmetry from the generator-level 
response. Small deviations in the MC closure are taken as a source of systematic uncertainty. 
Other systematic uncertainties have been estimated from the extrapolation procedure and from 
the scaling of the central-core widths between data and MC samples. The final results for the 
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scaling factors are presented in Fig. 40 These results demonstrate that, given the current data 
statistics, the observed data over MC ratios of the resolutions tails are within a factor of 1.5. 
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Figure 40: The data/MC scaling factors for the tails of the resolutions observed in the dijet 
samples for different rj and pT bins, using the > 2.5cr window. 

7.4.2 7 + Jet Measurement 

This method uses the pj balance in 7+jet events to estimate the non-Gaussian tails of the res- 
olution. It can be used separately for the low- and high-response tails because the response 
distribution is not symmetric as in the case of dijet measurement. The number of events out- 
side 2.5(7" range is counted in bins of Pj, to compare the resolution tails from the MC simula- 
tion to the level observed in data. An example distribution is shown in Fig. 41 (left) for the 
32 < pj < 52 GeV bin and central tj, and the ratio of the number of tail events in data and MC 
samples vs. pj is shown on the right, as function of pj. The available statistics in data do not 
allow to measure the resolution tails with precision, but a constant fit to the ratio is consistent 
w^ithin uncertainties w^ith the study using the dijet sample, and thereby provides a cross-check. 
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Figure 41: An example jet pj resolution function for the 32 < pj < 52 GeV bin (left); ratio of 
the number of tail events in data and MC samples vs. Pj (right). 
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A study of the jet energy response and the p^-resolution in the CMS detector has been pre- 
sented. The various measurements were performed using the 2010 dataset of proton-proton 
collisions at \/s = 7TeV corresponding to an integrated luminosity of 36 pb^ . Three differ- 
ent jet reconstruction methods have been examined: calorimeter jets, jet-plus-track jets, and 
particle-flow jets, clustered with the arA\-kj algorithm with a distance parameter R = 0.5. 

The jet energy response of all jet types is well understood and good agreement between data 
and simulation has been observed. The calibration is based on MC simulations, while residual 
corrections are needed to account for the small differences between data and simulation. The 
calibration chain also includes an offset correction, which removes the additional energy inside 
jets due to pile-up events. Various in situ measurements, which utilize the transverse mo- 
mentum balance, have been employed to constrain the systematic uncertainty of the jet energy 
scale. For all jet types, the total energy scale uncertainty is smaller than 3% for pj > 50 GeV 
in the region |?/| < 3.0. In the forward region 3.0 < \rj\ < 5.0, the energy scale uncertainty for 
calorimeter jets increases to 5% (Fig.[28|. 



The jet pr-resolution has been studied, using the dijet and 7+jets samples in both data and 
simulation. For PF jets in the region \rj\ < 0.5 with a pj of 100 GeV the measured resolution 
in the data is better than 10% (Figs. [34p5] |. The core as well as the tails of the jet pr-resolution 
function have been estimated, and close agreement is observed betw^een the 7+jets and dijet 
samples. The core of the measured jet pj-resolution in data is broader than the one obtained 
from the simulation, by 10% in the central region and up to 20% in the forward region. The 
resolution tails are in agreement with the simulation w^ithin statistical uncertainty. 
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